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On near Chandrasekhar mass central stars of planetary nebulae 
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ABSTRACT 

Aims. A recent spectroscopic analysis of central stars of planetary nebulae (CSPNe) claims that the sample includes five CSPNe with masses 
very close to the Chandrasekhar limit of white dwarfs. This claim should be verified or discarded from the available kinematical and chemical 
abundance information. 

Methods. Kinematical parameters are extracted from Galactic orbits and compared with parameters expected for populations of different ages. 
The chemistry of the nebulae is compared with average values for different types. 

Results. The reported high masses are not supported by our investigation. The claimed high central star masses are in contradiction with all 
other evidence. A more consistent picture emerges if CSPN masses close to the peak of the white dwarf mass distribution are assumed. 
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1. Introduction 

Stellar masses are the most fundamental parameter of stars. 
However, they are notoriously difficult to derive by direct meth- 
ods. The large majority of mass determinations make use of 
stellar models derived from computations of their evolution. 
Therefore it is of utmost importance that these models are 
checked against independent mass determinations. 

In this respect the situation for planetary nebulae (PNe) and 
their central stars (CSPNe) is mostly dire. Napiwotzki ( 1999) 
performed an NLTE model atmosphere analysis of high grav- 
ity CSPNe and determined masses from the post-AGB tracks 
of [Blocker (1995) and Schonberner (1983). In a subsequent 
study iNapiwotzkil (l200lh showed that the resulting distances 
are consistent with trigonometric parallax measurements and 
distances determined from companions in wide binaries, thus 
confirming the parameter and mass determination. However, 
this comparison was limited to high gravity central stars, which 
had akeady entered the white dwarf cooling sequence. No such 
test has been performed for low gravity CSPNe on the con- 
stant luminosity part of the post-AGB evolution. "Statistical" 
distance determinations using properties of the PN, like the fa- 
mous Shklovsky^nefliod, exist, but they are notoriously unreli- 
able (lNapiwotzkj|200 1|) and therefore cannot provide an inde- 
pendent test. 

In a recent studv lPauldrach et alJ ( |2004 hereinafter PHM ) 
performed a model atmosphere analysis of UV spectra of nine 
luminous CSPNe. This analysis is based on spherically sym- 
metric NLTE model atmospheres, which treat the hydrodynam- 
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ics of the stellar wind in a self-consistent way. These and sim- 
ilar model atmospheres have been successfully applied to the 
spectra of hot m assive main sequence stars (see discussion in 
IPauldrach et al.h2004) . The beauty of this method comes from 
the fact that it enables the determination of the absolute quan- 
tities stellar mass M and radius R without reference to stellar 
structure models. M and R are closely linked to the terminal 
wind velocity Voo and the mass loss rate M, which can be deter- 
mined from a fit of the spectral lines. 

The'PHM" analysis of the CSPNe sample produced a quite 
astonishing result (Table[0: five out of the nine analysed stars 
have a mass ver y close to the Chandrasek har limit for white 
dwarfs (l■40M,7^: lHamada & Salpeteill96ll) and a sixth one has 
a mass of 1 . IM©. This result is astonishing in at least two ways: 
1) white dwarfs of such a high mass are quite rare and thus we 
wouldn't expect such a high fraction in a sample of CSPNe, 
and 2) the masses derived by PHM differ quite considerably 
from what is derived from the comparison of effective temper- 
ature Teff and gravity g with theoretical post-AGB tracks. If 
the results of PHM can be proofed correct, this would indicate 
substantial flaws in today's theory of CSPNe evolution. 

Extraordinary claims need extraordinary proof. How can 
the PHM results be verified? One possible test is the compari- 
son of the derived spectroscopic distances with other indepen- 
dent measurements. However, as already discussed in IPHMI 
there is a deplorable lac k of an y reliable distance determi- 
nations for stars from the IphmI sample in particular and for 
CSPNe of similar type in general. 

In this article I will take a different pathway. The initial 
mass-final mass (IMFM) relation for white dwarfs is well es- 
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Table 1. Stellar parameters of the high mass CSPN from lPHMl heliocentric radial velocities (from iDurand et al.lll993) and 
proper motion measurements. The last column gives the source of the proper motion values. 



PNG 


common name 






M 


d 


Vrad 




yUa source 






(K) 


(cm s"^) 


(Me) 


(kpc) 


(kms-') 


(mas yr"') 


(mas yr"') 



025.3+40.8 IC4593 40000 3.80 1.11 3.63 22.0±0.5 -8.7±1.6 4.1±1.5 TY2 

083.5+12.7 NGC6826 44000 3.90 1.40 3.18 -6.2±0.6 -11.0+2.9 -9.7±1.7 TRC 

215.2-24.2 IC418 39000 3.70 1.33 2.00 61.9±0.5 -1.2±1.7 2.5±1.8 TRC 

315.1- 13.0 Hen2-131 33 000 3.10 1.39 5.62 -1.2±4.4 -2.1±1.6 -5.1±1.6 TY2 
316.1+08.4 Hen2-108 39000 3.70 1.33 6.76 -11.1±0.4 -0.6±2.7 -1.1±2.5 TY2 

345.2- 08.8 IC 1266, Tel 35 000 3.62 1.37 3.73 -84.1±4.7 -1.6±2.2 -11.4±2.2 TY2 



TY2: Tycho-2 catalogue 1H02 et al. 2000J; TRC: Tycho reference catalogue l H0g et alll99i 



tablished ('Weidemann"2000^. Some uncertainty exists at the 
high mass end, but we know from observations of white dwarfs 
in open clusters and binaries that high mass white dwarfs are 
produced by high mass progenitors. Even if one allows for 
some scatter in the IMFM relation and blue straggler scenar- 
ios it is hard to imagine that the progenitor of a > IM© CSPNe 
is older than 1 Gyr, only slightly shorter than the main sequence 
lifetime of a 2Mq star cSchaller et al..l992.) . Therefore the high 
mass CSPNe should obviously belong to the young thin disk 
population of the Milky Way. This population is characterised 
by small peculiar velocities (ie. a small velocity dispersion) 
and small scale heights perpendicular to the Galactic disk. 
Progenitors of high mass CSPNe are associated with dredge- 
up of elements produced by nuclear burning, which can be de- 
tected by an analysis of the CSPNe or the surrounding PNe. 

I will use the kinematical properties and the PN chemistry 
of the PHM sample for a sanity check. In Sect.|3]the kinematics 
of this sample is compared to the expectation for a young pop- 
ulation. Sect. 13 presents the evidence available from nebular 
abundance measurements. Sect.l^will conclude with a discus- 
sion of the combined evidence. 

2. Kinematics 

The kinematical state of a sample of stars contains informa- 
tion on the population membership and their age. The motion 
of stars in the Milky Way is usually described in the orthog- 
onal coordinate system X, Y, Z and the corresponding veloc- 
ities U, V, W. X points from the Sun in the direction of the 
Galactic centre, Y points in the direction of the galactic rota- 
tion at the position of the Sun and Z towards the north Galactic 
pole. Here I will assume 8 kpc distance of the Sun from the 
Galactic centre, an orbital velocity of the local standard of rest 
(LSR) of 220km s^'and Ue = 10.0km s"', Vq = 7.2km s-\ 
Wq = 7.2 km s^'relative to the LSR. 

Older populations are characterised by larger scale heights 
in Z direction and larger velocity dispersions in all directions. 
Investigations of solar neighbourhood stars make often use 
of diagrams plotting various combinations of the velocities 
U, V, W. However, this would produce misleading results, if ap- 
plied to stars several kpcs away from the Sun. 

Insights into the kinematical state of stars can be 
gained from their orbits in the Milky Way. The orbits 

S iresented here are calcu lated with the programme orbit6 
Odenkirchen & Broschd [T992) . The Galactic potential of 
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Fig. 1. Meridional cut of the Galactic orbit of IC418. The orbit 
is integrated over 5 Gyr. 



IC 4593 

ecc = 501 

= 5-5 kpc 
nze = 0.603 




Fig. 2. Meridional cut of the Galactic orbit of IC 4593. 



lAUen & Santiiiaiil (Il99ll) is adopted. Radial velocities and 
proper motions were compiled from litera ture a nd are pre- 
sented in Tabled Distances were taken from PHM'. The results 
for IC 41 8 and IC 4593 are shown in Figs.[2and|5] respectively. 
We use so-called meridional cuts with p as hori- 

zontal axis. Note that the orbit of a thin disk star without any 
peculiar velocities would appear as a dot in these plots. The or- 
bit of IC418 (Fig. ^ shows the characteristics of a disk star. 
The eccentricity ecc, defined as 

ecc^(R,-Rp)KR,+Rp) 
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Table 2. Adopted velocity dispersions cry, cry, ctw and asym- 
metric drift Vatj for the young thin di sk (YTD), the old thin disk 
(OTD) and for the thick disk (from iRobin et alJ2003l) . 



population 


(Tu 


(Tv 


CTw 


Vad 




(kms-') 


(kms-') 


(kms-') 


(km s-') 


YTD (0.15-lGyr) 


19.8 


12.8 


8.0 


3.1 


OTD(7-10Gyr) 


43.1 


27.8 


17.5 


14.8 


thick disk 


67 


51 


42 


53 



Fig. 3. Kinematical properties of the CSPNe in the eccentricity- 
nze diagram. The distances computed by PHM and measured 
proper motions from Table^were used for this plot. The con- 
tour lines show the 3cr limits for young disk, old disk and thick 
disk stars (in order of increasing size). The error bars for indi- 
vidual objects indicate 68% confidence intervals. 

with Rg, and being the apo- and perigalactic distances, is 
significant, but not extreme. The maximum distance from the 
Galactic plane Zmax is on the high side and makes member- 
ship in the old thin disk or the thick disk possible. The orbit of 
IC4593 (Fig.|2ji is more extreme reaching large distances from 
the Galactic plane and having high eccentricity. This points to 
a halo nature of this star. The discussion of bigger samples and 
error limits is facilitated by condensing the properties of the 
orbits into simple numbers. I will make use of the eccentricity 
ecc and normalised z-extent 

nze - Zmax/p(Zmax) 

introduced by Ide Boer et al.l lll997h to take the effect of the 
Galactic potential diminishing with galactrocentric distance 
into account. 

Values of ecc and nze for the CSPNe sample are plotted in 
Fig- El The measurement errors were propagated via a Monte 
Carlo simulation (distance errors of 15% were assumed). For a 
discussion of these results we need to know where to find pop- 
ulations of different age in this diagram. We started with the ve- 
locit y dispersions and asymmetric drifts given by Robi n et all 
( 1200 3) for the young thin disk, the old thin disk and the thick 
disk (Table |2j and translated this into distributions in ecc and 
nze via a Monte Carlo simulation with the orbit6 programme. 

The results are shown in Fig.|3l They are completely con- 
trary to what we expect for a young sample. Only one CSPNe 
(IC418) is marginally consistent with a membership in the old 




Fig. 4. Kinematical properties of the CSPNe in the eccentricity- 
nze diagram with proper motions treated as free parameters 
(filled squares). The parameters resulting from the observed 
proper motions are indicated by open squares connected to the 
new results. 

thin disk, Fig.|3lindicates even a halo nature of the supposedly 
young CSPNe of IC 4593 ! Admittedly, the error bars are large. 
The largest contribution comes from the uncertainties of the 
proper motion measurement. In principle the measurement of 
proper motions on this level of accuracy can be hampered by 
nebular structure close to the central star, which could cause 
offsets. If one wants to play devil's advocate one could argue 
that the proper motions are completely unreliable. Although 
this is very likely over-pessimistic, it is possible that the error 
bars underestimate the real errors. 

Thus I performed an experiment by treating the proper mo- 
tions as free parameters. Distance and radial velocities were 
kept fixed at their values from Table ^ and an optimisation 
was performed to minimise the combined values of nze and 
eccentricity, yeccP^+nzj? . The result is shown in Fig.|2 Even 
now only three CSPNe can be found close to the old thin disk 
contour, not to mention the young thin disk. Three stars (He 2- 
108, He 2-131, IC4593) must be members of the thick disk or 
halo populations with a very high level of confidence. Thus we 
conclude that the kinematical evidence is in clear contradiction 
with the CSPNe being part of a young population. 

3. Nebular abundances 

High mass central stars with high mass progenitors are usu- 
all y identifie d with PNe of type I in the classification scheme 
of IPeimbertI lll978l) . These show enhanced helium and ni- 
trogen abundances (He/H > 0.125 and N/O > 0.50; 
IPeimbert & Torres-Peimbert 1983!), which are explained by 
dredge-up processes in the progenitor Although it is by no 
means clear, whether all PN with high mass CSPN have en- 
hanced abundances of helium and nitrogen, evidence for a 

j 

type I nature of the lPHMl objects would strongly support their 
high mass nature. 

Chemical abundances determined from the analysis of the 
surrounding nebulae are compiled in Table|3l These were taken 
from the compilation of Perinotta C1991) unless noted other- 
wise. Table|3lis supplemented by average abundances of type I 
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Ta ble 3. N ebular abundances for the very high mass CSPN of 
the lPHMl sample. 



PN 


He/H 


N/0 


log(0/H) + 12 


IC4593 


0.099 


0.047 


8.50 


NGC6826 


0.107 


0.10 


8.46 


IC418 


0.086 


0.13 


8.76 


Hen2-131 




0.38 


8.67 


Hen 2-108" 


> 0.117 


0.28 


8.40 


IC 1266'' 


0.074 


0.09 


8.71 


type I 


0.137 


0.85 


8.66 


type II-III 


0.103 


0.26 


8.66 




The picture becomes very different, when we assume a low 
mass close to th e peak of the white dw a rf mass distribution 
(0.56 . . .O.59M0 'Naniw otzki et"Ii1ll999t iTJebert et al J 120051) 
and a radius consistent with standard stellar structure calcula- 
tions. To demonstrate the effect I comp uted distances resu lting 
from the O.565M0 post-AGB track of SchonbernerVl983') and 
repeated the orbit calculations with these values. Note that I 
could have adopted the parameters of iKudritzki et alJ (Il997h 
instead, but it is not impossible that these analyses have there 
own caveats, as discussed by the authors. The result is shown in 
Fig- El Now the parameters of five CSPNe are consistent with 
being members of the old thin disk. IC 4593 maybe member of 
the thick disk or halo. Interestingly, the oxygen abundances of 
this PN is somewhat lower than normal (Tabled), which would 
be consistent with this finding. On the other hand the kinemat- 
ical properties of single objects are no firm proof of their affil- 
iation. They could have, for example, suffered from particular 
violent encounters with other stars. 

We conclude that neither the kinematical evidence nor the 
PNe chemistry s upport s the claim of high masses of the six 
CSPNe from the'PHM' sample. The findings are far better ex- 
plained, if one assumes much lower CSPNe masses, which 
would correspond to lower progenitor masses. Unfortunately 
this would indicate that the model atmospheres applied by 
|PHM are not yet able to produce reliable results for CSPNe, 
in contrast to their successful application to massive stars. 



Fig. 5. Kinematical properties of the CSPNe in the eccentricity- 
nze diagram resulting for assumed CSPNe masses of 0.565Mq 
as discussed in the text. 
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and type II-III PNe from IPerinottol lll99lh . PNe of type II 
and III are identified with lower mass CSPNe from the Galactic 
disc, probably including some objects from the thick disk as 

well. 

I Q] l i'x that none of the PNe of the lPHMl sample 
qualifies as type I. The chemical abundances are well within 
the range expected for run of the mill PNe of type II and III. 
We conclude that again a high mass nature of the investigated 
objects is not supported. 

4. Putting it all together: What is the true nature of 
these objects? 

This article presented two tests to verify or disprove the results 
of|PHM that six of the CSPNe analysed by them have very high 
masses. The kinematical results of Sect. |2 are clearly at odds 
with these CSPN being young objects with massive progenitor 
stars. The PNe chemistry pr esented in Sect. |3l also points to 
lower mass progenitors of the lPHMl CSPNe. Viewed together I 
conclude that the PHM sample has failed this test and the high 
masses are not confirmed. 

Now, are all these stars members of the old and metal poor 
thick disk and halo populations? This seems unlikely, given that 
the PN abundance determinations in Table |3l show that these 
are not very metal poor objects. However, the halo-like kine- 
matics of the,WlM sample are largely the result of the very 
large distances computed from the PHM parameters (Tabled). 
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